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Abstract 

We demonstrate that the detection of shock modulations of the neutrino spectra 
from a galactic core-collapse supernova is sufficient to obtain a high significance de- 
termination of the neutrino mass hierarchy if the supernova event is observed in both 
a Mton-class water Cherenkov detector and a 100 kton-class liquid argon detector. 
Neither detailed supernova neutrino flux modelling nor observation of Earth matter 
effects is needed for this determination. As a corollary, a nonzero value of 6 X will be 
established. 
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1 Introduction 



The current status of neutrino oscillation parameter estimations can be very briefly sum- 
marized PP as follows: Atmospheric (solar) neutrinos oscillate with \8m 2 a \ ~ 0.002 eV 2 and 
9 a ~ 7r/4 {Sm 2 s ~ 8 x 10" 5 eV 2 , 9 S ~ tt/6 jSj) 1 . All we presently know about 6> x . is that 
sin 2 9 X < 0.05 at the 2a C. L. jlj. A long-standing hope is that neutrinos from a core-collapse 
supernova (SN) may shed light on two of the unknown oscillation parameters, sgn(<5m 2 ) and 

Only a handful of neutrinos from a Type II SN have ever been detected. The detection 
of 11 neutrinos from SN 1987A in Kamiokande II [5] and 8 neutrinos in the Irvine Michigan 
Brookhaven experiment fH] have been of great importance for understanding core-collapse 7J. 
It is evident that the physics potential offered by a future galactic SN event is immense. 
With cognizance of this potential, experiments dedicated to SN neutrino detection have 
been proposed jH] even though only a few galactic SN are expected per century. 

Attempts have been made to extract neutrino oscillation parameters from the 19 SN 
1987A events. However, conclusions drawn from such analyses are highly dependent on 
the neutrino flux model adopted and are far from robust. For example (and within the 
context of this paper), it was claimed that the data favor the normal hierarchy (5m 2 . > 0) 
over the inverted hierarchy (5m 2 < 0) provided sin 2 6 X > 10~ 4 0, but this conclusion was 
contradicted in Ref. jTUj. 

Neutrinos from a galactic SN could in principle provide a wealth of information on neu- 
trino oscillations. A determination of 9 X and the neutrino mass hierarchy from SN neutrinos 
is unique in that ambiguities fTJ ^] arising from the unknown CP phase 5 and the devia- 
tion of atmospheric neutrino mixing from maximality do not corrupt it. The absence of the 
eight-fold parameter degeneracies that are inherent in long baseline experiments J2| results 
because (i) nonelectron neutrino fluxes 2 do not depend on 5 independently of neutrino con- 
version [T31, an d so SN neutrinos directly probe 6 X , and (ii) whether atmospheric mixing is 
our notation, 5m 2 s (Sm 2 ) is the solar (atmospheric) mass-squared difference and 8 S , 9 X and 9 a are the 

mixing angles conventionally denoted by 9\2, 613 and 623, respectively £Q. 

2 We focus on detection via charged current u e and v e interactions, which cannot distinguish between the 

different nonelectron neutrino species (that we denote by v x with x = fj,, r, ft, f ). 
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maximal or not is immaterial since 8 a does not affect the oscillation dynamics. 

Investigations of the effect of neutrino oscillations on SN neutrinos in the context of a 
static density profile (i.e., neglecting shock effects) have been made in Refs. [TH H5j . 
Whether or not the mass hierarchy can be determined and 9 X be constrained depends sensi- 
tively on the strength of the hierarchy between (E Pe ) and (E Vx ). The higher (E Ux ) / (Ep e ) 
is above unity, the better the possible determinations [14j. Unfortunately, modern SN 
models that include all relevant neutrino interaction effects like nuclear recoil and nucleon 
bremsstrahlung indicate that the hierarchy of average energies is likely smaller than ex- 
pected from traditional predictions; (E Ux ) / (E Pe ) is expected to be about 1.1, and no larger 
than 1.2 as opposed to ratios above 1.5 from older SN codes jTHj. Another relevant 
uncertainty is that different SN models predict different degrees to which equipartitioning 
of energy between z/ e , v e and v x is violated. For example, in Ref. ^H] an almost perfect 
equipartitioning is obtained while according to Refs. [201 121] equipartitioning holds only to 
within a factor of 2. 

Given these uncertainties, it is not a simple task to determine 9 X and the mass hierarchy 
simultaneously from SN data (16 . A significant improvement would be a determination of 
the mass hierarchy independently of predictions for (E Ux )/(E Pe ) and equipartitioning from 
SN models. In this paper we propose a new method using SN neutrinos to determine the 
mass hierarchy that exploits recent advances in the understanding of shock propagation in 
SN. 

At densities ~ 10 3 g/cm 3 , neutrino oscillations are governed by 5m 2 and sin 2 9 X |14j . 
Neutrinos (antineutrinos) pass through a resonance if 8m 2 a > (5m 2 , < 0). As the shock tra- 
verses the resonance, adiabaticity is severely affected causing oscillations to be temporarily 
suppressed, as first pointed out in Ref. [221- After the shock moves beyond the resonance, 
oscillations are restored. Then one expects a dip in the time evolution of the average neutrino 
energy and the number of events 3 . This modulation is visible in the neutrino (antineutrino) 
3 Recently, this idea was taken one step further in Ref. |23]. A reverse shock caused by the collision 
between a neutrino-driven baryonic wind and the more slowly moving primary ejecta may also form. The 
direct and reverse shocks yield a "double dip" signature |28j . In the present work we restrict our attention 
to the effects of the forward shock which is a generic feature of SN models and whose existence is better 
established than that of the reverse shock. 
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channel for a normal (inverted) mass hierarchy and only if tan 2 9 X ^> 10~ 5 i.e., only for oscil- 
lations that would occur adiabatically for a static density profile. Previous work exploiting 
the dip to obtain information about oscillation parameters can be found in Ref. 124*1 • 

Within the first 4 seconds or so, the violation of adiabaticity caused by the shock is felt 
only by neutrinos with energy less than about 20 MeV. At these energies some models predict 
the v x flux to be larger than the v e and v e fluxes (2nj and others predict the converse jTUj . 
At later times, the shock affects higher energy neutrinos for which all models predict the v x 
flux to be dominant. The dip is observable even for {E Vx ) = {E Ve ) because the fluxes are 
flavor-dependent (23]. Thus, a signature in high energy neutrinos a few seconds after bounce 
is quite model-independent. Through our analysis we show that the signal is so robust that 
a restriction to high-energy events is unnecessary. 

We investigate the significance with which the mass hierarchy can be determined by 
measurement of the u e spectrum at a 100 kton liquid argon detector and the v e spectrum at 
a 1 Mton water Cherenkov detector from a galactic SN at a distance of 10 kpc with binding 
energy 3 x 10 53 ergs. The number of unoscillated events in the liquid Ar (water Cherenkov) 
detector is expected to be O(10 5 ) (O(10 6 )). (Although a liquid argon detector can measure 
both the v e and v e spectra, it is an order of magnitude more sensitive to the v e flux than to 
the v e flux [23]). We are interested in the detectability of a dip in the time evolution of either 
the v e or the v e spectrum, but not both. By correlating the two spectra, it should be possible 
to establish the mass hierarchy and that tan 2 9 X ^> 10 -5 . The presence of a dip will also 
provide further evidence that SN simulations correctly depict shock propagation. If a dip is 
not found in either channel, it would suggest that tan 2 X < 10 -5 , since the salient features 
of the theory of core-collapse SN have already received strong support from SN 1987A jjj. 

2 Shock density profile and neutrino oscillation prob- 
abilities 

Realistic time-dependent density profiles of SN are obtained from detailed numerical siu- 
mulations. As one moves in towards the neutron star, the profile at a given instant has a 
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sharp density rise followed by a rarefaction region where the density can drop significantly 
below that at the outer edge of the shock. The authors of Ref. [22] have provided a generic 
time-dependent density profile that mimics those of supernova simulations. The shock front 
is steepened artificially to reintroduce the physical requirement of a density discontinuity 
which is lost in hydrodynamic simulations due to the limited (few 100 km) spatial resolu- 
tion. We adopt the empirical parameterization of this profile (which is continuous in the 
supernova radius and time) developed in Ref. [2o] . 

Because of the dip in density in the rarefaction region, neutrinos may hop between mass 
eigenstates up to 3 times before leaving the SN envelope. Under the assumptions that 
the transitions factorize and the neutrino phases can be averaged away, a simple analytic 
expression for the overall hopping probability Pff(|5m„|, 9 X ) was obtained in Ref. [2E], which 
agrees remarkably well with phase-averaged results of Runge-Kutta evolution of the neutrino 
flavor propagation equations. We employ the analytic expression for Pjj to calculate the u e 
and v e survival probabilities. These survival probabilities do not depend on Sm 2 since we 
are neglecting Earth matter effects [Hj: 

sin 2 S cos 2 6 X P H + sin 2 6 X {1 - P H ) , (1) 

cos 2 6 S cos 2 9 X , (2) 

sin 2 9 S cos 2 9 X , (3) 

cos 2 6 S cos 2 6 X P H + sin 2 X (1 - P H ) . (4) 

Here, the N and I subscripts denote normal and inverted hierarchy, respectively. We note 
that the factorization of the 3-neutrino dynamics into two 2-neutrino subsystems continues 
to hold with the subsystem governed by Sm 2 and 9 S remaining adiabatic (for the now well- 
established Large Mixing Angle solution |3]) as in the case of a static density profile j2H]. 



P N (u e -> u e ) = 

Pn(v £ -> u e ) = 

Pj{y e -> v e ) = 

Plipe -> Ve) = 



3 Neutrino spectra 

We use the parameterization of Ref. for the primary unoscillated neutrino spectra given 
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where i = u e , u e , u x , and (Ej) (/%) is the average energy (dimensionless shape parameter 
that quantifies the width of the spectrum) of species i. In principle, both (Ei) and /3j can 
be time- dependent. Throughout, we assume f3 Ve = 4, (3p e = 5 and (3 Vx = 4 fTTj. Here, 
is the neutrino emission rate or luminosity of species i. The only available simulation that 
tracks neutrino emission for times long enough to enable studies of shock effects on neutrino 
oscillations is that of the Livermore group. We adopt the luminosities of Ref. [TU] . 

Our goal is to demonstrate the model-independence of our new method. To this end, we 
consider models that span a wide range of predictions for the energy spectra. We conserva- 
tively consider the following parameter ranges for the initial spectra (indicated by the "0" 
superscript) in our analyses (with all energies in MeV): 

14<«><22, «) = (0.5-0.9)«), «> = (1.0-1.6)«>, (6) 

0.5 < $° e /< < 2.0 , 0.5 < $° e /< < 1.6 . (7) 

This is in line with the hierarchy of average energies (E® ) < (E®) < (E® ), expected from 
the well-known interaction strengths of neutrinos with matter; the species which interacts 
less decouples earlier and has a higher temperature. 

We will occasionally refer to a specific model which appears to be quite pessimistic for 
determining oscillation parameters because the different neutrino species have very small 
spectral differences. It is a model obtained by the Garching group by accounting for all 
relevant neutrino interaction effects. The initial spectra predicted by the Garching model 
have [20] 

«> = «> = 15 MeV ' < E ° c > = 12MeV > </< = 0-5, </< = 0.5. (8) 

4 Simulations 

For the 1 Mton Hyper-Kamiokande detector j27j we assume an energy threshold of 7 MeV 
and a fiducial volume for SN neutrinos of 685 kton which is consistent with the fiducial to 
total volume ratio expected for the UNO detector jS]. We only consider events from inverse 
/3-decay and use the higher-order cross section of Ref. [28J. We assume that the energy 
resolution is that of the Super-Kamiokande detector 
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For the 100 kton liquid Ar detector with near-perfect efficiency, we consider the 
absorption process, v e + 40 Ar — > e~ + 40 K*, whose cross section is calculated in Ref. 
The effective threshold is 6 MeV. Since the energy resolution is expected to be extremely 
good, we do not include resolution effects. 

For each mass hierarchy and several values of 9 X , we randomly choose SN models defined 
by the parameter ranges of Eqs. (|b17|) and simulate 10, 000 possible v e and v e spectra at the 
two detectors (including statistical fluctuations). We fix tan 2 6 S = 0.4 [3] and \5m 2 a \ = 0.002 
eV 2 .2]. We ignore Earth matter effects (thereby not committing ourselves to specific zenith 
angles for the SN) which only become comparable to shock effects at very small 9 X . 

We illustrate the effect of the shock on the neutrino spectra in Fig. ^ 

5 Statistical analyses 

We consider two cases (a) (Ei) = (E®) and (b) with (Ei) falling linearly with time 5 sees 
post-bounce according to (E^t > 5)) = (£°)(1. 04167 - 0.008333t). We have chosen this 
parametrization to agree with that in Ref. j2H|. Although this fall-off in energy with time 
is arbitrary, it serves to demonstrate the robustness of our analysis to specific assumptions 
about (Ei(t)). After all, in a real SN, nature will make a unique selection for (Ei(t)) which 
is as yet unknown to us. 

So long as tan 2 6 X 3> 10 -5 , we expect a dip in the v e (u e ) spectrum if the hierarchy 
is normal (inverted) as in Fig. ^ A way to establish the channel in which the dip occurs 
exploits the fact that the dip is expected to occur 4 seconds after the shock forms (i.e. after 
the core rebounds). We split each spectrum into 2 time bins. All the events occuring in 
the first 4 seconds constitute the "early time bin" , with mean energy (E etb ) and the events 
occuring between 4 and 10 seconds constitute the "late time bin", with mean energy Em,. 
Their ratio is 

R= (E ltb )/(E etb ). (9) 

In the case that (Ei) = (Ef), for the normal hierarchy we expect R(u e ,Ar) to be smaller 
than unity and R(v e , H^O) to be close to unity and vice-versa for the inverted hierarchy. If 
tan 2 X < 10~ 5 , we expect both ratios to be close to unity. 
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For falling (Ei), we expect both R(v e , Ar) and R(v e , H 2 0) to be slightly below unity even 
if tan 2 9 X < 10~ 5 . This is because (Ei) starts falling 5 seconds after bounce, which lowers 
(Eitb) relative to (E et b). This shock-independent suppression cannot be large since even at 
t = 10 sees, (Ei) = 0.96(Ef). (Clearly the size of the effect depends on how fast (Ei) falls). 
Then, the hierarchy is deduced by comparing the relative deviations of R from unity. If 
R(p e ,H 2 0) < R(v e ,Ar) < 1, the hierarchy is inverted. If R(v e , Ar) < R(H e ,H 2 0) < 1, the 
hierarchy is normal. 

In Fig. |2K, we plot R(v e , Ar) vs. R{v e , H 2 0) for the normal hierarchy (red squares) and 
for the inverted hierarchy (green triangles) with tan 2 9 X = 0.01 assuming (Ei) = (Ef). The 
blue dots are for 9 X = 0. There is little overlap between the clusters, and the inverted 
hierarchy can be established at more than 3cr. As can be seen in Fig. |2b, the case in which 
(Ei) fall linearly also yields little overlap between the clusters, although the clusters are now 
systematically displaced to lower values of R(v e , Ar) and R(v e , H 2 0). 

We now focus on the (Ei) = (E®) case. It is no surprise that it is more difficult to 
identify the normal hierarchy considering the lower statistics in the Ar detector. However, it 
is possible to establish the normal hierarchy reasonably well. In Fig.Et the dashed red curve 
shows the fraction of spectra corresponding to the normal hierarchy and tan 2 6 X = 0.01 with 
R(v e) Ar) below a given value. For example, 88% (98.7%) of the spectra corresponding to 
the normal hierarchy and tan 2 9 X = 0.01 have R(v e , Ar) < 0.98 (0.99). The solid blue curve 
shows the number of spectra with 9 X = and R(v e , Ar) below a given value divided by the 
number for the normal hierarchy and tan 2 9 X = 0.01 with R(v e , Ar) below the same value. 
Thus from Fig. Et, the fraction of spectra with 9 X = that mimic the spectra involving 
shocks, and have R(v e , Ar) < 0.98 (0.99), is 0.33% (7.6%). Said differently, 88% (98.7%) of 
the spectra indicate the normal hierarchy correctly with a probability of 99.67% (92.4%). 

We emphasize that these results have included the Garching model of Eq. (JHJ). As ex- 
pected, the results of a similar analysis of 1, 000 spectra simulated from the Garching model 
alone are somewhat worse. Specifically, we find that 87% (99%) of the spectra simulated 
with the normal hierarchy and tan 2 9 X = 0.01 have R(v e , Ar) < 0.98 (0.99), and identify 
the normal hierarchy with a confidence of 99% (87%). It is remarkable that even for initial 
spectra with tiny differences, such as in the Garching model, it is possible to establish the 
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normal hierarchy. 

Figure |3b is similar to Ot, but for the inverted hierarchy. Clearly, a determination of the 
inverted hierarchy is easier. We see that 99.5% of the spectra simulated with an inverted 
hierarchy and tan 2 9 X = 0.01 have R(u e , H2O) < 0.99, and indicate the inverted hierarchy at 
the 99.99% C. L. 

The case in which (Ei) fall linearly with time 5 seconds after bounce gives similar results 
as those for the constant (Ei) case. 

We have also investigated if the hierarchy can be determined for tan 2 9 X = 0.001. It is 
still possible to establish the the inverted hierarchy with high-confidence. For the normal 
hierarchy we find that while shock effects remain visible, it is no longer possible to establish 
the normal hierarchy with high-significance. 

6 Conclusions 

Our analysis has established that detections of shock modulations in either the u e spectrum in 
a kton-class liquid Ar detector or in the v e spectrum in a Mton-class water Cherenkov detector 
will provide definitive evidence for the true neutrino mass hierarchy and that tan 2 9 X 3> 10 -5 . 
We emphasize that the modulations should be observed in one and only one of the two 
channels. 

A similar analysis can be performed supposing that a reverse shock also develops in 
the SN envelope (221 ■ Then, one expects to see a double-dip signature in either the u e 
or v e spectrum if 9 X is not too small J2H1- in addition to probing oscillation parameters, a 
presence of a double dip will provide confirmation that SN simulations correctly depict shock 
propagation. 

We have supposed that both a large water Cherenkov and large liquid argon detector are 
available. A possible alternative to a large water Cherenkov detector is the IceCube detector 
since good energy resolution is not a requirement of our analysis. Then R(P e , Icecube) can 
be compared with R(u e ,Ar) 

An aspect that bears further consideration is whether our analysis is robust under vari- 
ations of the shock density profile. We are currently investigating this issue in Ref. [32J. 
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Figure 1: The expected v e (solid green curves) and v e spectra (dashed red curves) for the 
normal hierarchy and inverted hierarchy, respectively with tan 2 X = 0.01. The thick (thin) 
lines correspond to (E^) = (Ef) {{E,j) that fall linearly 5 sees after the core-bounce). The 
curves are based on the (pessimistic) Garching model of Eq.(JEJ). Note the pronounced dips 
in the v e (z/ e ) spectra for the normal (inverted) hierarchy. 
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Figure 2: (a) R{y e ,Ar) vs. i?(V e , H%0) for SN models spanned by Eqs. (jtiJ7j) and tan 2 ^ = 
0.01 assuming; (Ei) = (Ef). R is the ratio of the mean energy of the events occuring between 
4 and 10 seconds to the mean energy of the events occuring in the first 4 seconds. The 
red squares (green triangles) correspond to the normal (inverted) hierarchy. The blue dots 
are for 9 X = 0. Each contour contains 95% of the 10 5 v e and v e spectra generated for 
the corresponding case. We have plotted only 10 3 triangles, squares and dots so as to not 
overwhelm the figure, (b) is the analog of (a) but for the case in which (Ei) falls linearly 
with time 5 sees after bounce. The separation of the clusters is a measure of the significance 
with which the hierarchy can be determined. 
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Figure 3: (a) The dashed red curve shows the fraction of spectra corresponding to the normal 
hierarchy and tan 2 9 X = 0.01 with R{v e) Ar) below a given value (use the y-scale on the left). 
The solid blue curve shows the number of spectra with 9 X = and R{v e) Ar) below a given 
value divided by the number for the normal hierarchy and tan 2 9 X = 0.01 with i?(z/ e , Ar) 
below the same value (use the y-scale on the right), (b) is the analog of (a) but for the 
inverted hierarchy. Spectra with (Ei) = (E®) were assumed. See Eq. for the definition 
of R. 
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